Quantum pump effect in one-dimensional systems of Dirac fermions 
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We investigate the behavior of the directed current in one-dimensional systems of Dirac fermions 
driven by local periodic potentials in the forward as well in backscattering channels. We treat the 
problem with Keldysh non-equilibrium Green's function formalism. We present the exact solution 
for the case of an infinite wire and show that in this case the dc current vanishes identically. We also 
investigate a confined system consistent in an annular arrangement coupled to a particle reservoir. 
We present a perturbative treatment that allows for the analytical expressions of the dc current in 
the lowest order of the amplitudes of the potential. We also present results obtained from the exact 
numerical solution of the problem. 
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I. INTRODUCTION 

The possibility of employing ac-fields to induce di- 
rected flow of charge and spin in quantum systems is 
in the limelight of quantum transport phenomena. Re- 
cent experimental observations of the quantum pumping 
effect, includingthe periodic deformation of the walls of 
quantum dots [U, 0, [1, Q and the induction of currents 
by applying surface acoustic waves in carbon nanotubes 
[5| triggered an important development in the theoret- 
ical descrmtion of time-dependen t q uantum trans por t 
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The development of new technologies based on quan- 
tum phenomena, is also gauging the search of novel ma- 
terials to achieve transport in devices of reduced scale. 
Among the many structures under investigation, devices 
fabricated of graphite are capturing an increasing atten- 
tion. From the theoretical side, the Physics of graphite 
has the appealing feature of being a realization of a gas of 
Dirac fermions in reduced dimensions. In fact, the behav- 
ior of the conductance of single- wall nanotubes [1^ seems 
to be reasonably explained on the basis of a Liittinger 
liquid, which is one of the few models of interacting elec- 
trons in condensed matter Physics that is amenable to 
be analytically solved. The Physics of graphene, being a 
single graphite sheet seems also to be suitably described 
as a two-dimensional gas of Dirac fermions [24|. These 
developments are motivating a renewed interest in the 
Physics of Dirac fermions and in the application of quan- 
tum field theoretical methods to shed light on low di- 
mensional Condensed Matter problems. In particular, 
a precise understanding of the effect of time-dependent 
perturbations on these systems will be crucial in order to 
build nanodevices based on the pumping of charge and 
spin into a conductor. Recently, we have solved the prob- 
lem of an ac-potential in a one-dimension infinite wire of 
Dirac fermions, emphasizing the separate role of the for- 



ward and backward channels in the resulting total energy 
density and transport properties [25| . 

The breaking of certain symmetries has been singled 
out as a necessary condition to generate directed cur- 
rents. This idea has been firstly introduced in the con- 
text of classical "ratchet" systems, [26'] and also seems 
to apply in paradigmatic models of quantum pumps, like 
two barrier structures driven by two harmonic fields os- 
cillating with a phase-lag. [2l| We recall that the pump 
effect is characterized by the induction of a directed cur- 
rent with pure ac driving, i.e. without introducing any 
explicit dc bias voltage. In this work we adapt the two 
oscillating barriers setup to a one-dimensional wire of 
Dirac fermions. One of our aims is to identify which are 
the minimal ingredients to achieve the quantum pump ef- 
fect in these systems. An infinite wire of Dirac fermions 
with two local forward ac fields oscillating with a phase 
lag can be analytically solved by generalizing the treat- 
ment of Ref. [2a As we will show, the ensuing induced 
current has a vanishing dc component. Then, it is nat- 
ural to explore the role of confinement in relation to the 
appearance of non vanishing dc currents. This is the 
main purpose of this work. In particular, we consider a 
finite wire of length L bent into a ring and coupled to an 
external reservoir while driven by ac forward and back- 
ward fields, as indicated in the sketch of Fig. [1] This 
setup enables us to analyze the role of dynamical for- 
ward and backscattering in wires containing a spectrum 
with many resonant levels, being in contact with a dissi- 
pative environment. We treat the problem with Keldysh 
non-equilibrium Green's functions [2^ and find analyti- 
cal expressions for the dc current at the lowest order of 
perturbation theory in the pumping amplitudes. We also 
present a system of exact equations and exact results for 
the dc current obtained by the numerical solution of these 
equations. 

The paper is organized as follows. In sections II and III 
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FIG. 1: (Color online) Sketch of the setup. Above: Infinite 
wire driven by periodic fields inducing dynamic forward and 
back scattering at the positions x — ±a/2. Bellow: The ac 
fields act on a finite wire of length L in anular geometry. The 
system is coupled to a particle reservoir a the position x = 0. 



we present the model and the theoretical treatment. In 
section IV we present the exact evaluation of the dc cur- 
rent in an infinite wire with dynamical fowrd scattering. 
Section V is devoted to study the transport behavior of a 
ring with finite length coupled to a reservoir. A summary 
and conclusions are presented in Section VI. Technical 
details, in particular, the procedure followed to obtain 
an effective action through the functional integration of 
the reservoir's fields, are described in Appendix A. 



II. MODELS 



We study a one-dimensional system of Dirac fermions 
in two different setups. Firstly, we consider an infinite 
wire with two harmonically time-dependent local poten- 
tials applied at the positions xj = ±a/2 (see Fig. la). 
This system is described by an action 



S — >S'o + Sac 



(1) 



where Sq is the linearized free dispersion relation for right 
and left-movers (in our unit system h = vp = e = 



So = ^J dxdt{4{t){dt+dMt)+i^}it)idt-d,)Mt)), 

(2) 



Sac describes forward and backward scattering processes 
through the two time-dependent potentials: 



Sac = dxdtSix-Xj)\vfit)i4{t)ijrit) + 



(3) 



Note that, for simplicity, the dependence of the fields 
on X has been omitted in the above equations. To be 
specific we consider harmonic potentials with equal am- 
plitudes and a phase-lag 6: V^'^ (t) = V^'^ cosftot and 
V.f'^{t) = V^'^ cos {Hot + S), where F, B denote forward 
and backward channels, respectively. 

In the second configuration, we assume that the Dirac 
fermions propagate in a ring of length L, which is in con- 
tact to a fermionic reservoir. Also in this case, harmon- 
ically time-dependent local potentials are applied at the 
positions Xi = (L±a)/2 . The action for the full system 
can be casted in terms of four terms: 



S — So + Sac + Sr 



Sc 



(4) 



The reservoir is assumed to be an infinite system of free 
fermions which couples to electrons in the ring through 
a forward scattering term. The term Scant account for 
the interaction of fermions with the reservoir, placed at 
a; = 0. It reads: 

Scout = j dxdtvc6{x){i^l{t)xr{t) + i^\{t)xi{t) + 

+ximr{t) + xlmm, (5) 

where the fields xj(i), denote degrees of freedom of 
the reservoir. The latter can be integrated out, giving 
rise to an effective self-energy So(t — t') (see Appendix 
A) which leads to the following effective action for the 
internal degrees of freedom of the ring: 



Seff = j dx dt dt' (jpl{t) 

iidt + d^)5{t - t') ~ Six) So(i - i')] Ait') 
+i;j{t) [i(dt - d.,)S{t - t') - 5{x) Eo(< - <')] Mt')) 

2 

-|- / dxdt5{x — Xi) 

2—1 

^[vf{t){i,l{t)Mt) + 4>\{t)Ut)) 
+v,''{t){^l{t)Mt) + i^\{t)Mt))\ (6) 

Our main purpose is to compute the dc component of 
the time-dependent charge current. The latter is defined 
through the condition of the charge conservation: 



dpi{x,t) dpr{x,t) di{x,t) 



dt 



dt 



dx 



(7) 
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being pj{x,t) = {ipl^{x,t)'ip^{x,t)) , j = l,r, which casts: 
j{x,t) = pr{x,t) - pi{x,t). (8) 
The corresponding dc current is 

■ rlr r 



r"- = — I dtj{x,t), 

To Jo 



(9) 



where Tq = 27r/r2o is the period of the osciUations. Note 
that, due to the conservation of the charge, (see eq. (O), 
it does not depend on the spacial coordinate x. 



III. THEORETICAL TREATMENT 

We solve the problem presented in the previous section 
by recourse to Keldysh non equilibrium Green's func- 
tions. This formalism is based on the solution of the 
equation of motion of the matricial Green's function ob- 
tained from the corresponding action. 

From actions (HJ and ^ it is easy to show that the 
matricial Green's hmction Q reads; 



G++{x,x',t,t') G+-,{x,x',t,t') 
G-+(x,x',t,t') G--,{x,x\t,t') 



where the indexes -j-, — refer to the closed time path of 
Keldysh formalism [23, being 

G++,{x,x' ,t,t') = -i{T^^{x,t)i;\,{x\t')), 

G--ix,x',t,t') = ~i(f^^ix,t)il;\,ix',t')), 

G+-(x,x',t,i') = t{^l,{x',t')i^^{x,t)), 

G-+{x,x\t,t') - ~t{ij^ix,t)i^l,ix',t')), (11) 

where T and T denote temporal and anti-temporal or- 
dering, respectively. These functions satisfy the following 
Dyson equation: 

(^[i{dt + Sjdx) - V^{x,t)]Sj.^^ - 

-V^{x, t)Sy^^j-^^gj-^^^>{x, x', t, t') 

—5{x)5^^^^ j dtiY:o{x,t ~ ti)Q^^^j'{x,x' ,ti,t') 
= 5^,^>5{x-x')5c{t-t') (12) 

where Sc{t~t') is the matricial delta function (extended 
to the time contour): 

6{t~t') t,t' on G+ 
Scit-t')={ -S{t~t') t,t' on G- (13) 
otherwise 

V^'^ix,t) ^S{x-xi)Vf'^{t)+6{x-X2)Vi''^it), = 



In terms of Green's functions, the current ([8]) reads 

j{x,t) = -i[GX'{x,x,t,t) - Gji^{x,x,t,t)]. (15) 

As shown in the next section, for the infinite system 
with So(i — t') — 0, the Green's function admits a fac- 
torization in terms of the noninteracting (equilibrium) 
Green's function, at least for ac potentials that behave 
as forward scatterers. One then obtains an analytical 
expression for j{x,t) with vanishing time average. This 
result can be extended perturbatively to the backscatter- 
ing case. 

The situation changes dramatically when fermions 
propagate in a ring connected to a reservoir. As a con- 
sequence of a non-instantaneous term associated to the 
coupling to the reservoir, the equation of motion for the 
Green's function ceases to be exactly solvable even for 
the forward channel and its solution must be found by 
recourse to perturbation theory or numerical methods. 

IV. INFINITE WIRE 

For the infinite wire, we can adapt the treatment of 
Ref. [1^, to the case of more than one oscillating po- 
tentials. For the case of pure forward ac potentials (i.e. 
Vj^ — 0,V«), j{x,t) can be exactly evaluated and the re- 
sult is a purely ac current. To show this we solve the 
Dyson equation (|12p for 'So{t — t') = and two oscillat- 
ing barriers placed at x = ±a/2. Specializing eq. (jl2p 
to the case = and considering first, for illustrative 
purposes the Green's function for right-movers, results 

it{dt + d.,) - V^{x, t))g{x, x',t, t') = 5{x - x') 5c{t - t') 

(16) 

It is very useful to note that the solution of the above 
differential equation for the matrix Q can be obtained in 
terms of the free function Q'^ , the latter being the solution 
of (HH), with V^{x,t) = 0. The result is: 

g{x, t- x', t') = g°{x- x', t - t') exp[/3(a;, t) - P{x' , t')], 

(17) 

with 



/3(x,t) 



dx' dt' G^'^ix -x',t- t'^V'^ix', t'), (18) 



where G^{x -x',t~ t') = G°'++ - G°'+- is the equilib- 
rium retarded Green's function of the free system. The 
above equation can be readily solved to obtain: 

I3{x,t) = -iV^{6{x-a/2) cos 17o(i - (a: - a/2)) + 

+e{x + a/2) cos {na{t - (x + a/2)) + (5)Vl9) 



for the right and left movers, respectively, and 

s++(t-t') s+-(i-0 

Eo(<-t')=| I- (14) 

So +(*-*') So-(^-<') 



Having now /3(x, t) we have, via equation (jl7p . every com- 
ponent oiQ. Of course, the same steps can be followed to 
evaluate the Green's functions for left-movers. The func- 
tions G^~ {x, a;', t, t') can be used to compute the current 



4 



from eq. (fT5|) . This leads to 

jix,t) = ^^^(0ix + a/2)sin[Qo{t~ix + a/2))] 
Zn \ 

+ 9{x - a/2) sin[17o(t - (x - a/2)) + S] 
+ e{-{x + a/2))sm[no{t + x + a/2)] 

+ 9{-{x - a/2)) sm[no(t + x- a/2) + (5])(20) 

which yields = 0. 

In the case 7^ it is not possible to arrive at a 
closed expression for the current, one is forced to perform 
a perturbative expansion. We have computed j(x,t) up 
to second order in , showing that even for this case 
the time average of j{x, t) vanishes for the infinite wire. 



V. FINITE RING IN CONTACT TO A 
RESERVOIR 

A. Evaluation of the Green's function 

In the case of the annular setup, the presence of a self- 
energy in the effective action S^f / renders the procedure 
of the previous section unsuitable even in the case of 
pure forward dynamical scattering. However, the Green's 
functions can be formally exactly calculated by following 
a similar route to the one adopted in the solution of time- 
dependent lattice Hamiltonians [13] ■ In what follows we 
summarize this procedure. It is convenient to combine 
the different components of the Greens function matrix 
to define a retarded Green's function: 



G^^^^.{x,x\t,t') = Q{t-t')[G-+,{x,x\t,t') 



(21) 



and to perform a Fourier transform with respect to the 
difference of times: 



G ^ ^/ (x, X ^ t^t ) 



(22) 

Substituting it in the integral representation of the 
Dyson's equations ([1^. the following set of linear equa- 
tions is obtained [2(j |: 

G^^y{x,x\t,Lo)^G"^ix,x',uj) 
2 

+ ^ e-'^''°*-'^^G!^^y {x, xj ,t,u; + no) 



Gy{xj,x',uj) 



+ e(*"«*+^^)G^,y (x, X, ,t,u- no) 



G'l,{xj,x',uj), 



(23) 



being r = I and I = r. The function G^[xj,x' ,w) is the 
retarded Green's function of Dirac left and right movers 



of the ring in contact to the reservoir with Vi{t) = 0. 
This function can be calculated by considering the corre- 
sponding Dyson equation, which in the present case takes 
the following simple form: 

G°{x,x',uj) = g°{x,x',uj) 

+G°(x,0,c^)E«HgO(0,x',c^), (24) 

where g^{x,x\u!) are the free retarded functions, corre- 
sponding to the ring without driving and uncoupled from 
the reservoir, while T,q{uj) is the Fourier transform of the 
retarded self-energy: E^{t-t') = e{t - t' )[!:„+ {t - t') - 
T,Q~ (t — t')]. It can be in general expressed in terms of its 
spectral representation r(cj) = — 2Im[E^(ijj)] as follows: 



duj r('i^) 



2tt uj — uj + iO+ 
The equation for the lesser function reads: 

G+-,{x,x',t,t') = J2 

ll=l-X 



(25) 



xS+-(c.)[G«^^(0,a;',i',c.)]*. (26) 



This function can be used to evaluate current from p5|) . 
Recalling that 'E,'^~{uj) — if{uj)T{ui), where /(w) — 
9{—{u} — ji)) is the Fermi function (we shall consider the 
zero-temperature case, /x is the chemical potential of the 
reservoir), the dc current can, thus, be expressed as fol- 
lows: 



•dc 



duj 
2^ 



fi^mco), 



(27) 



where the functions T{u!) collects the spectral contribu- 
tions to the current from the forward and backward chan- 
nels. Concretely: 



T{lu) = - / " dt(\ Gf,(x,0,t,oj) r 
To Ja V 



G§{x,0,t,iu) 



G«(x,0,i,c^) 



G«(x,0,i,c.) I r(c.) 



(28) 



Therefore, the evaluation of the retarded Greens func- 
tions leads us to the current. To this end, we must solve 
Dyson's equation 



B. Perturbative expansion in the pumping 
amplitudes 

In view of the obvious difficulties of solving equation 
we propose in this section a perturbative solution 
that allows us to evaluate to order 0{\V^\'^, \V^\'^). 
We will obtain an analytical expression for j'^'^, through 
the "transmission function" T(lu), valid for any kind 
of fermionic reservoir interacting via forward scattering 
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with the fermions in the ring, hke the one described by 
Scant introduced in eq. ([5]). 

Before going to the perturbative treatment, let us note 
that given the Green's function of the free system: 

, , n expiuj{x - x' - L/2) 
(X, x,u)^g,{x,x,u.)= \un^ ^—^^ 

(29) 

Eq. pi]) can be exactly solved: 



G'^ix, X ,uj) ^ G\'{x ,x,uj) ^ hm — , 

A-»oo A hlyUj) 

(30) 



with 



hiuj) ^ 2 sincjL/2 - iE^{Lo))e-"^^/^ . (31) 



The normalization factor tt/A is introduced because the 
functions G'^{x,x,uj) and g^{x,x' ,uj) are not regular, 
in the sense that the densities of states pj{x,x,uj) = 
— 2Im[G!^(x, x, w)] and /9° (x, X, cj) — —2lm[g^{x,x,u!)] do 
not integrate to 1. Actually, the ensuing integrals diverge 
and, for this reason, upper and lower cutoffs ±A are in- 
troduced to ensure the correct sum rule. In the case of 
the Green's function for the ring coupled to the reservoir, 
the normalization is the same as for the free system pro- 
vided that the reservoir has a finite bandwidth W, such 
that T{w) = for |w| > W, and A » W. 

Given a concrete model for the reservoir, which im- 
plies assuming a concrete functional form for the spec- 
tral density T(uj), one has a definite function h{uj) (See 
([3T|l ). This, in turn, leads to the retarded functions for 
the system without time-dependent potentials, given in 
([M)) . These functions are our building blocks for the 
perturbative solution of the full retarded Greens func- 
tion GJ^^{x,x' ,t,uj). Solving equation ([^ up to second 
order in the couplings and replacing the results in 
we get 



iS X 



where 

rW(a 



A2 I hico) I 

sin noa[T^^\Lo)+T^^\Lo)] 

-yJ[r(3)(c.) + rW(^)]y 



1 



h(LU + ilo) 



h{uj~ no) p' 



(32) 



(33) 



T'-^\uj) = 2Re 



-zOn-L/2 



T'^^\u;) = 2Re^—( 
^ ' h*(u;)\ 



h* (lo) \h* {lo + no) /i*(tj-f7o) 



e*"-^ /sin {{2uj - no)a)eJ^"^^^ 



(34) 



(llj) \ h*{iLj + no) 

sin((2w-f f7o)a)e-''^^'"^/2 
h*{uj - Qq) 
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FIG. 2: (Color online) The function T{uj) for pure forward 
scattering. Parameters are: Vf ~ 0.2, 5 = n/2, a — 2, 
L = 40. We consider a reservoir with bandwidth VK = 16 
and a coupling constant Vc — 0.1. Upper and lower panels 
correspond to fio ~ 0.1,0.15, respectively. 



r(4)(a 



sin (2a; — f2o)a sin(2ci; + rio)a 



h{uj — no) 



hiio -f no) 



(35) 



In equation (|32p we have omitted, for sake of simplicity, 
the limit A ^ oo. Let us stress again that these formu- 
las do not depend on the specific choice of the function 
r(a;). This allows us to obtain some general properties 
of quantum pumping in the annular geometry which are 
independent of the details of the reservoir. As a first 
point, we stress that unlike the case of the infinite wire, 
in the finite ring, forward scattering is enough to generate 
a directed current. Second, it is interesting to note that 
we recover the behavior j'^'^ oc sin(J) predicted in Ref. 

for a two-barrier quantum pump in a linear setup, 
also obtained within Floquet scattering matrix formal- 
ism [l^ as well as within Green's function formalism 
in tight-binding models [13, [111- As a third point, for 
(5 ^ modiji) and fioa ^ 1, the dc current behaves like 
j'^'^ cx r^Oj which is the typical behavior of the so called 
"adiabatic" regime of quantum pumps [3, [HI, [IBi [l3i [3 ■ 
A final and important point to note is the role of geomet- 
rical factors, like the length L and the situation of the 
barriers a. In particular, the prefactor in (j32p indicates 
that for weak amplitudes and special combinations of pa- 
rameters such that noa = tt, the dc current is inhibited 
in the forward channel. In addition, geometrical factors 
lead to sign modulation of the induced dc current as a 
function of the chemical potential /i, as we will discuss in 
more detail in the next section. 
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FIG. 3: (Color online) dc current J'*'^ as a function of the 
chemical potential of the reservoir for pure forward scattering. 
Parameters are the same as in Fig. [2l 
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FIG. 4: (Color online) The function T{u!) for pure backward 
scattering with Vb = 0.2. Other details are the same as in 
Fig. El 



C. Exact numerical solution 

In this section we present results obtained from the 
numerical exact solution of the Dyson's equation 
In order to avoid the problem of the regularization of 
the Green's function discussed in the previous section, 
we consider in this section functions 

n— — N 

being A;„ = 2tt/L while the upper (lower) sign corre- 
sponds to 7 = r, L This function coincides with (|29p in 
the limit of iV — > oo, and has well defined spectral prop- 
erties, i.e. it satisfies the usual sum rule for the spectral 
density. The drawback of this Green's function is that 
Dyson's equation ([24| cannot be analytically solved. In 
any case, the equilibrium Green's functions can be nu- 
merically evaluated from the expression: 

G'°(x,x',cj) = g°(x,a;',Lj) 

+ 1-E«(c.)g0(0,0,c.)' ^ ^ 

and used as input to the numerical evaluation of (j23[) . 

So far, we have not introduced any specific model for 
the reservoir. In the numerical procedure we have repre- 
sented the latter system in terms of a semicircular density 
of states of bandwidth W and a coupling Vc between ring 
and reservoir. The ensuing self-energy is defined from 
(Uni) with T{uo) = A\vc\'^Q{W - \uj\)^/W^ - /W'^ . 

Results for the function T(w) in the case of pure for- 
ward scattering are shown in Fig. [2l This function dis- 
plays a landscape of peaks and anti-peaks that are in 
accordance with the structure predicted by the pertur- 
bative solution. For Vb = only the functions T^^^((jj) 



and r(2)(a;) contribute (see Eqs. jMl) and JM])). 

Notice that these functions have peaks centered at the 
zeros of the function h[ijj ± VLq) with So(w) = 0, which 
correspond to the positions of the energy levels of the un- 
coupled ring, i.e. at uj±^l^ = 27171 / L, with n integer. The 
modulation with period ujL of T^'^^{uj) also explains the 
pattern of alternative maxima and minima of T{lj) at the 
energy levels of the uncoupled ring. In the lower panel of 
Fig. [2] it is shown the behavior of T(w) for fip ^ 27r/L, 
i.e. close to resonance with the characteristic frequency 
associated to the level spacing of the uncoupled ring. 

At zero temperature, the integral of T{lo) over energies 
below the chemical potential gives the dc component of 
current for a given jjL. The behavior of J'''^ as a function 
/i is shown in Fig. [3] for a driven system with the same 
parameters of [51 In the upper panel, corresponding to 
VLq < 27r/L, a sequence of plateaus of width rjp are clearly 
distinguished. In addition sign inversion of the current 
takes place at points distant in the level spacing 2-k/L. 
In the case of the lower panel (close to resonance) the 
plateaus approximately coincide with the level spacing 

The function T{uj) for pure dynamical backscattering 
{Vp — 0) is shown in Fig. [4l The behavior of this func- 
tion is also in agreement with the one predicted by the 
perturbative solution. In this case, only the functions 
T^^\lij) and T'-^\uj) contribute. As in the case of pure 
forward scattering, these functions have a sequence of 
peaks and anti-peaks defined from uj±D,o = 27r/L. How- 
ever in the present case, the structure is richer than in 
the case of forward scattering. This is due to the exis- 
tence of two modulating factors: e*"^ in T^^\lo), and the 
sm[{2Lj±no)a\ in T(3)(w) as well as in T('*)(tj). While the 
first factor introduces a periodicity which is coincident 
with the level spacing, the second one introduces a peri- 
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FIG. 5: (Color online) dc current J'*'^ as a function of the 
chemical potential of the reservoir for pure backward scatter- 
ing. Other details are the same as in Fig. U 



odicity na that is defined by the position at which the ac 
potentials are applied. In addition, for backscattering, 
different terms contribute with different signs, while in 
the case of forward scattering all the contributions are 
additive. 

In Fig. m an envelope with period air can be clearly 
distinguished and it is indicated with an arrow in the 
lower panel. The corresponding dc current is shown in 
Fig. [5l For low pumping frequencies (upper panel) as 
well as for frequencies close to resonance (lower panel), 
the behavior of the directed current as a function of the 
chemical potential displays sign inversions and a sequence 
of maxima and minima separated in 7ra/2. In addition, 
there is a fine structure of plateaus of width fio, which 
is more pronounced in the case of low flo (upper panel). 
Also notice, and additional structure related to the level 
spacing energy 27r/L of the ring. 

Thus, dynamical forward and backscattering mecha- 
nisms leave clearly different fingerprints in the behavior 
of the induced directed current as a function of the chem- 
ical potential of the reservoir. While both processes are 
sensitive and give rise to patterns related to the level 
spacing and the pumping frequency, backscattering leads 
to a response that also provides information on the posi- 
tions at which the pumping potentials are applied. This 
feature accounts for interference processes in the elec- 
tronic propagation that take place when electrons emit 
or absorb quanta and invert at the same time their di- 
rection of motion at the backscattering centers. 



VI. SUMMARY AND CONCLUSIONS 

In this work we have studied quantum pumping in- 
duced by two local ac potentials in the forward and back- 
ward channels, oscillating with a single frequency and a 



phase lag, in a one-dimensional system of Dirac fermions. 
This problem can be exactly solved analytically in the 
case of an infinite wire and pure forward scattering. The 
exact solution shows that the charge current in this case 
is purely ac without a dc component. For this reason, we 
have investigated a confined system of annular shape in 
contact to a particle reservoir. We have employed a for- 
malism based on Keldysh non-equilibrium Green's func- 
tions to solve this problem. In this case, an analytical 
solution is available a the lowest order in the pumping 
amplitudes, while for arbitrary amplitudes, it can be ex- 
actly solved by numerically evaluating a Dyson equation 
for the retarded Green's functions. 

We have shown that a directed current establishes in 
the confined geometry with dynamical forward as well 
as backscattering. As a function of the phase-lag 6, the 
current shows the cx sin (5 behavior already presented in 
the literature in the framework of adiabatic descriptions 
3, models of plane- wave ^16] and tight-binding electrons, 
20, HH and experimentally observed in Ref. [J. Interest- 
ingly, in the case of pure forward scattering there is an ad- 
ditional factor oc sin(f7oQ)j denoting interference effects 
related to the geometrical arrangement. As a function of 
the chemical potential of the reservoir, forward and back- 
ward scattering display features in the behavior of the dc 
current which are related to the wire level spacing the 
pumping frequency. In the case of backscattering, there 
is additional structure related to geometrical parameters. 
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APPENDIX A: INTEGRATION OF THE 
RESERVOIR'S DEGREES OF FREEDOM 

In this appendix we show how to get the effective action 
given in eq. ^ . To this end we consider the path-integral 
representation for the piece of the partition function that 
involves the fields Xr and xi to be integrated: 



dxl dxl dXr dxi exipi{Sres + Scont) (Al) 



where 



Sres = J dxdt {xlDr Xr 



+ xlDrXi) (A2) 



where and Di are differential operators which specific 
forms depend on the dispersion relations satisfied by the 
particles in the reservoir. Since the integration proce- 
dure can be performed quite generally, we will specify 
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the actual form for these operators at the end of the 
derivation. On the other hand, the interaction between 
the reservoir's modes and fermions in the ring is given by 
Scant (see eq.(IS|)). In order to obtain an action which is 
quadratic in the reservoir's fields we perform a transla- 
tion in these fields of the form Xr.i ~* (^r.i such that: 



Of course, there are also two terms depending quadrat- 
ically on ipr,i which will also contribute to the effective 
action. The general form of these action terms is: 



Xr.i (a;, t) = arj (x, t) + o^,; (x, t) , 



(A3) 



(and similar transformations for xl i)- Choosing the func- 
tions ar.i as 



arj{x,t) — —Vc J dx' dt' D^i{x,x';t,t') S{x')ipr.i{x' ,t'), 

(A4) 

(and similar expressions for a^^^) the resulting action is 
quadratic in the ar,i fields, which are then integrated out. 



dx dt dt' iPl^i {x, t) 5{x) D;! (0, 0; t, t') tpl^i {x, t). 

(A5) 

In this paper we will consider particle reservoirs such 
that right and left propagators coincide for x = x' . We 
will also assume time translation invariance such that 
D;^{{),Q]t,t') = D-^{Q,{)]t - t'). The reservoir's self- 
energy used in the paper is defined in terms of this prop- 
agator as Y^o{t - t') = vlD^^{Q,Q;t - t'). An explicit 
form for this function is chosen in section IV C. 
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